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Abstract — A discrete-time two-user interference cliannel model 
is developed that captures non-linear phenomena that arise 
in optical liber communication employing wavelength-division 
multiplexing (WDM). The effect of non-linearity is that an 
amplitude variation on one carrier induces a phase variation 
on the other carrier. Moreover, the model captures the effect of 
group velocity mismatch that introduces memory in the channel. 
It is shown that both users can achieve the maximum pre-log 
factor of 1 simultaneously by using an interference focusing 
technique introduced in an earlier work. 

I. Introduction 

The increase in traffic demand and the advances in optical 
technology over the past two decades made determining the 
capacity of optical fiber networks of great interest. The optical 
fiber channel suffers from three main impairments of different 
nature: noise, dispersion, and Kerr non-linearity. The interac- 
tion between these three phenomena makes the problem of 
estimating the capacity challenging. One approach is to solve 
the propagation equation numerically and then approximate 
the capacity by estimating the input and output statistics (see, 
e-g-j CI and references therein). A second approach is to 
estimate the fiber capacity analytically, but then one must make 
simplifying assumptions. In 12J, we studied a simplified non- 
linear model in which all orders of dispersion were neglected. 
In this paper, we add first-order dispersion to the non-linear 
model, which makes the model more realistic. 

II. Fiber Channel Models 

Suppose that two optical fields at different carrier frequen- 
cies uji and Ld2 are launched at the same location and propagate 
simultaneously inside the fiber. The fields interact with each 
other through the Kerr effect |f3] Ch. 7]. Specifically, neglecting 
fiber losses, the propagation is governed by the coupled non- 
linear Schrodinger (NLS) equations ||3j p. 264, 274]: 
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where i = \/— T, Ak{z, T) is the time-retarded, slowly varying 
component of field k, k — 1,2, the f32k are group velocity dis- 
persion (GVD) coefficients, the 7^ are non-linear parameters. 



and d = /3i2 — /3ii where the (3ik are reciprocals of group 
velocities. The parameter d is called group velocity mismatch 
(GVM). To simplify the model, we assume that the second 
order dispersion effects are negligible, i.e., /32i — /322 = 0. 
In previous work ||2l, the group velocity mismatch was also 
ignored (d — 0), which leads to a memoryless channel model. 
In this paper, we consider the case of non-zero GVM {d ^ 0), 
which introduces memory as we shall see. The coupled NLS 
equations ([l]l-(|2]) have the exact solutions |3 p. 275] 



A2(L,T) ^ A2{0,T ^ Ld)e">'^^^'^'> , 
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where L is the fiber length, and the time-dependent non-linear 
phase shifts are obtained from 



ML,T)^-fi[\A,{0,T)\^L + 2J \A2{0,T-zd)\^dzy 
ML, T) = 72 {\A2{0, T)\^L + 2 / |Ai(0, T + zd)\''dz^ . 



where z = is the point at which both fields are launched. 
Kerr non-linearity does not change the shape of the pulse in 
the time domain but causes interference through intensity- 
dependent phase shifts, which affects the spectrum of the 
pulse. The non-linear phase shift seen by a field due to 
the field itself is called self-phase modulation (SPM), i.e., 
7fci|^fe(0, r)|2 is an SPM term, whereas the non-linear phase 
shift seen by a field due to other fields, which causes interfer- 
ence, is called cross-phase modulation (XPM). 

In ||2|, a discrete-time model was obtained by sampling 
the fields at the transmitters and the receivers. In general, 
this introduces noise with large variance at the receivers 
because wide-band noise is not removed by filtering. In the 
next section, we overcome this shortcoming by using matched 
filtering before sampUng at the receivers. 

III. Continuous-Time Model 

Consider the case of d 7^ 0. Without loss of generality, 
suppose that d > 0. Let {xk['m]}m=o ^^ '^^e data sequence 
sent by transmitter k. Suppose that the transmitters employ 



square pulse shaping, i.e., the signal sent by transmitter k is 



^fe(0,T)=^XfcHp(r-mr3). 



where 



,^. ^ r ves/Ts, < i < t, 

^ ' \ 0, otherwise. 

The signal observed by receiver k is 

rkiT)^Ak{L,T) + Zk{T), k = 1,2. 
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We write 0i(L, r) as 

0i(L,T) = 0n(i,T) + ,^i2(L,r), (16) 

in which we have defined 

(Aii(L,i)=7ii|^i(0,i)l', (17) 

^,2iL,t) = 2j,Lw — |A2(0,A)pdA, (18) 

^s Jt-Ld 

where L^f = Ts/d. Lw is called the walk-off distance. 
Since p(A) = for A ^ [0, T^], we have 



where Zk{T) is Gaussian white noise with E(zfc(r)) = 0, and 
E(zfe(r)z*(T + r)) == iV(5(T). The processes zi{T) and Z2(r) 
are independent. 

The analysis of the two receivers is similar, hence we 
focus on receiver 1. The received signal is fed to a bank 
of linear time-invariant (LTI) filters with impulse responses 
{hf{T)}f<,r^, where 
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/i/(i)-p*(-t)exp(i2^/VT,), 
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where ■(/'ig (0 i^ defined as 



and Fi C Z = {. . . , —1, 0, 1, . . .}. In other words, the bank 
of filters is a set of frequency-shifted matched filters whose 
impulse responses are orthogonal. The choice of the set Fi is 
specified in section [V] 

The output of the filter with index / is 

2/i,/(T)=ri(T)*/i/(T), 

where * denotes convolution. The "noiseless part" of the 
output of the filter with index /, yij{T), is given by 



47\t)^^f \p{\-mT,)\'dX. (21) 

Assume that Ld = MT^ for some M e {1, 2,3,.. .}. The 
larger the group velocity mismatch is, the larger M is, the 
more memory the system has. For Ld = MTg, we obtain 



4™^(t) = f|x< 



yij(T) = (ri(T)-zi(T))*/i/(T) 
= A,{L,T)*hf{T) 
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t — mTs, rriTs <t< (to + \)Ts 

Ts, (to + l)r, < t < (to + M)Ts 

(to, + M + l)Ts - t, 

(to + M)Ts <t<{m + M + l)Ts 
0, otherwise. 

(22) 



= J2 -^iH (p{T - TOT,)e'*^(^'^)) * hfiT) 



For r e [jTs^jTs + T,], we hav^ 
0ii(L,r) ^jiL^\xi[j]f, 



(23) 



E *^ 
(12) <Pi2{L,t) = 2-/iLw^[Y.\x2[j - r]\^ + 



= ^ xi [to] / p{t - mTs)p* (r - T) 

m=0 •^-°° 

e^ML,r)^-^2.f(r~T)/n^^_ (13) 

By sampling the output signal yij{T) at the time instants 
T = jTs, j e {0, 1,2,.. .}, we have 

yijiJTs) = X! ^1 ['"] / P^'^ ~ mTs)p* (t - jTs) 



(k2b-]p-ix2b--M]n^^^ 



(24) 



Hence, the integral in ( fTS) evaluates to 



I T, e^, ifB = 0, 
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where 
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Since p(t) = ^/EJT's for i e [0, T^] and p(i) = for 
t ^ [0, Ts], equation ( 14 1 reduces to 



'T, 



A = ^/^n|.Tl[j]|2 + z/ti2^ |a;2[j - r]\^ , 

S = ihi2{\x2[j]? - |a;2[j - Min - 27r/ 
/ill = liL^, and /ii2 = 27iLw— ^. 
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We use the convention of setting tlie quantities that involve a negative time 
index to zero. 



Therefore, the "noiseless part" of the output of the fiher 
with index / at time j, yi.f[j], is given by 



M 



mjb] = 

Xi[j]Es exp \ihii\xi[j]\'^ + ihu'^ \x2[j - 

r=l 

where 

r exp (z27r(7;i[j] -/))-! 
uij[j] = l i2Tr{vi[j] - f) 

in which we have defined vi[j] as 



(29) 



uijiJl, 



complex Gaussian random variables with mean and variance 
TV. We regard the h^i as channel coefficients that are time 
invariant and known globally. The following power constraint 
is imposed: 



n 



(38) 



j=i 



otherwise. 



V,[3\=hi2{\x2[3\?-\X2[3-M]\'')l2^. 
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The output of the filter with index / at time j, yij[i], is 

2/1, /[j] = yijiJTs) = yijlj] + zijij], (32) 

where 2:1./ [j] = zi{T) * hf{T)\T=jTs- The variable 21./ [j] is 
Gaussian with mean and variance NEg. Moreoever, due to 
the orthogonality of the filter bank impluse responses, we have 
E(5;i,/[j]z* f[j]) — for all f ^ f, which implies, because 
of Gaussianity, that the random variables {^i./[J]}/g.Fi are 
independent. 

IV. Discrete-Time Model 

Transmitter k sends a string XJ^ = 

{Xk[l],Xk[2],- ■ ■ ,Xk[n]) of symbols while receiver k 
observes Y^ = (r,[l], r,[2], • • • , F Jn]). The input Xk[j] 
of transmitter k to the channel at time j is a scalar, whereas 
the channel output l^i,[j] at receiver k at time j is a vector 
whose components are yfe,/[j], / G -Fk- The input-output 
relations are: 

Y,j[j]=X,[j] exp{thn\X,[j]\^ + 

M 

i/11251 l^2[j - r]f)Uij[3] + Zi.fij], (33) 

r=l 

M 

Y2j[j] = X2[j] exp (ih2i J2 \Mj + 2M - r]f+ 



A scheme is a coUection {{Ci{Pi, P2, N),C2{Pi, P2, N))} 
of pairs of codes such that at (Pi, P2, N), user k uses the code 
Ck{Pi,P2,N) that satisfies the power constraint and achieves 
an information rate Rk{Pi,P2,N) where fc = 1,2. We make 
a distiction betweeen two limiting cases: 1) fixed noise with 
growing powers, and 2) fixed powers with vanishing noise. 

Definition 1: The high-power pre-log pair (^1,^2) is 
achieved by a scheme {{Ci{Pi, P2, N),C2[Pi, P2, N))) if the 
rates satisfy 

fJ/3,7V)= lim Rk{PuP2,N) f j,^^ 2, (39) 

^ ' P.^oo:P.=Pf logiPk/N) 

where /3 is real. 

Definition 2: The low-noise pre-log pair (^11^2) is 
achieved by a scheme {{Ci{Pi, P2, N),C2{Pi, P2, N))} if the 
rates satisfy 



,(Pi,P2)= Im 



Pfe(Pl,P2,iV) 



for fc== 1,2. 



(40) 



N^Q logiPk/N) 

The (high-power or low-noise) pre-log pair (1/2. 1/2) can 
be achieved if both users use only amplitude modulation. In 
this work, we show that the high-power pre-log pair (1,1) can 
be achieved for any real /3 and positive N through interference 
focusing, which was introduced in |2|. Determining the maxi- 
mal achievable low-noise pre-log pair is a subject of ongoing 
work. 

V. Interference Focusing 

We use interferencing focusing, i.e., focus the phase inter- 
ference on one point by imposing the following constraints on 
the transmitted symbols: 



h2i\Xi[j]\^ = 2nni, ni = 1,2,3, 
hi2\X2[j]\^ = 2TTfi2, n2 = 1,2,3, 



(41) 
(42) 



^h22\X2[j + M]\'\U2j[j] + Z2j[j], 



(34) 



where 

Ukj[j] 



exp it2Tr{Vk[j] -/))-! 
i27riVk[j] - f) 



otherwise. 



(35) 



in which we have defined 14 [j] as 

Vi[j] ^ h,2{\X2[j]\' - 1X2 [j - M]\^)/27r, (36) 

^2[j] = h2ii\Xi[j + 2Af]|2 - |X2[j + A/]n/27r. (37) 

Zk.f[j] models the noise at receiver k at time j, the random 
variables {Zkj[j]}k,f.j are independent circularly-symmetric 



which ensures that the XPM interference is eliminated. 

Because of other constraints, e.g., power constraint, only a 
subset of the allowed rings is actually used, i.e., we choose 

1^1 b1 1' e T'l = {27rni//i2i, ni G Mi} (43) 

|^2b1|' e ^2 = {2Trn2/hi2, n2 & AA2} (44) 

where Wi, AA2 C N = {1, 2, 3, . . .}. In this case, Vk[j] e Vk, 
for k = 1,2, where 



Vk 



{mi - TO2 : mi G Wa^fe, m2 G Afs-k} ■. 



(45) 



which leads us to choose the sets of "normalized frequencies" 
Ti and J^2 of the filter banks at the receivers as J^i = Vi and 

-^2 = V2. 



Thus, under interference focusing, the output at receiver 1 
is a vector Y^^, whose components are {i^i,/}/eVi> whera^ 



Fij = Xi exp (ihnlXif) Uij + Z^j, 



and where 



Uij^ 



1, if Vi = /, 
0, otherwise. 



(46) 



(47) 



This means that, under interference focusing, exactly one 
filter output among all the filters "contains the signal" cor- 
rupted by noise, while all other filters put out noise. 

We wish to compute the information rate Ri — I{Xi]Y_^). 
From the chain rule, and non-negativity of mutual information, 
we have 



i?i=/(Xi;yi)>/(Xi;Yi,yJ 



Z 



l.Vi 



(48) 
(49) 



It was shown in ||2J that using interference focusing, we 
have 



lim 



>1, 



(50) 



Pl^oo log(Pl) 

which implies that ri > 1. 

Therefore, the (high-power) pre-log f i = 1 is achievable. 
Similarly, one can show that f2 = 1 is also achievable, and 
hence, the pre-log pair (1,1) is achievable. Using a simple 
genie-aided argumenr] we find that the maximal pre-log pair 
is (1,1). 

The following toy example illustrates our receiver structure, 
and the role interference focusing plays in choosing its param- 
eters. 

Example: Suppose that hi2 = 5, /121 = A, Pi — 8, and 
P2 = 7. Suppose that the users choose the sets of rings Mi — 
{1, 4, 9} and J\f2 — {2, 8} (see Fig.fTli, i.e., the power levels are 
Vi = {0.57r, 27r, 4.57r} and V2 = {0.87r, 3.27r}. These choices 
satisfy the power contraints and eliminate the interference. The 
parameters of the filter bank are J^i = Vi = {—6,0,6} and 
7^2 = ^2 = {— 8, — 5, — 3, 0, 3, 5,8}. In other words, receiver 

1 has 3 filters whose frequency responses are sine functions 
centered at /i — 6/Ts, /i, and /i + 6/Ts, whereas receiver 

2 has 7 filters whose frequency responses are sine functions 
centered at 7 different frequencies (see Fig. [2]). This shows 
that, because of the non-linearity, the receivers need to extract 
information from a "bandwdith" larger than the "bandwidth" 
of the transmitted signal. 

VI. Conclusion 

We introduced a discrete-time two-user interference chan- 
nel model based on a simplified optical fiber model. We 
assumed that second order dispersion is negligible. However, 
we considered non-zero group velocity mismatch as well 
as non-linearity. We showed that our discrete-time model is 
justified by using square pulse shaping at the transmitters and 

^We drop the time index for notational simplicity. 

'The genie gives each receiver the message it is not interested in. 




Fig. 1. Ring modulation used by transmitter 1 (left) and transmitter 2 (right). 
The thin lines are the rings allowed by interference focusing, and the thick 
blue lines are the rings selected for transmission. 




/2 - 8/r, 



Fig. 2. Frequency responses of the filters at receivers 1 (top) and 2 (bottom). 



a bank of frequency-shifted matched filters at the receivers. We 
proved that both users can achieve a high-power pre-log of 1 
simultaneously by using interference focusing, thus exploiting 
all the available amplitude and phase degrees of freedom. 
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